Structure, Structural Relaxation and Ion Diffusion in Sodium Disilicate Melts 
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We have investigated Na2Si2 0s melts with inelastic neutron scattering at temperatures up to 
1600 K. The Si-O network relaxes on a time scale of ns, whereas the Na ion relaxation dynamics 
are found on a time scale of 10 ps. The elastic structure factor exhibits at ~0.9 A -1 an emerging 
prepeak which becomes more pronounced with increasing temperature. The prepeak is caused by 
the formation of sodium rich regions in the partially disrupted Si-O tetrahedral network. 
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Silicate melts have great relevance both in earth sci- 
ence and technology: The physical properties of magma 
(molten silicate rock) dominate many geological pro- 
cesses and technological glasses are synthesized from the 
molten state. These silicate melts, natural and techno- 
logical, are multicomponent systems. A considerable ef- 
fort has been made to investigate a wide range of sil- 
icate melts with the general aim of linking structural 
with physical properties and to develop an atomic level 
understanding of their structure and dynamics [Q . Nev- 
ertheless, experiments on the structure and microscopic 
dynamics in the viscous melt well above the conven- 
tional glass transition temperature T g are rare. Here, 
we present inelastic neutron scattering results on sodium 
disilicate melts and show that at temperatures well above 
Tg significant changes in the elastic structure factor occur 
that are caused by sodium rich regions in the partially 
disrupted Si-O tetrahedral network. The existence of 
sodium rich regions has a strong influence on the sodium 
ion transport and may very well be the cause for the weak 
temperature dependence of the viscosity at higher tem- 
peratures. Our results on structure and dynamics can be 
rationalized using results of recent molecular dynamics 
simulations 

Alkali silicates are a simplified analog for most mag- 
mas and technological glasses. Compared with pure sil- 
ica, with a caloric glass transition temperature T g at 
~ 1500 K, the addition of Na20 partially disrupts the 
S1O2 network structure resulting in a halving of the T g 
at sodium disilicate. In addition, sodium disilicate has 
been the subject of extensive molecular dynamics simula- 
tions j2[jj] and thoroughly investigations HH- Inelastic 
neutron scattering covers a dynamic range that gives in- 
sight not only into dynamics on microscopic time scales 
but also into the medium range structure of silicate melts. 
In addition, the results can provide input for simulations 
used e.g. to model magma properties in volcanic activ- 
ity and are an experimental test of molecular dynamics 
simulations. 

Sodium disilicate Na2Si205 was synthesized from ul- 



trapure Na2C03 and Si02 powders. The mixture was 
melted in a Pt crucible and stirred with a PtsoRh^o spin- 
dle at superliquidus temperatures until it was homoge- 
neous and bubble free. Differential scanning calorime- 
try yielded a glass transition temperature at ~ 741 K in 
agreement with the literature value at T g ~ 737 K 0. 
For the neutron scattering experiment the sample was 
encapsulated in a Pt sample cell giving an annular sam- 
ple geometry of 40 mm in height, 22.5 mm in diameter 
and a 1.25 mm wall thickness. 

We performed inelastic neutron scattering measure- 
ments on the time-of- flight spectrometer IN 6 at the In- 
stitut Laue-Langevin in Grenoble. An incident neutron 
wavelength of A = 5. 1 A yielded an energy resolution of 
SE = 93 /ieV (FWHM) and an accessible wave number 
range at zero energy transfer of q = 0.3 — 2.0 A . Spec- 
tra were measured in the glass at 400 K and in the viscous 
melt at 800 K and between 1200 K and 1600 K in steps 
of 100 K. During cooling measurements were repeated at 
each temperature. No change between spectra obtained 
during heating and cooling could be detected. Additional 
high-energy resolution spectra at 300 K, 1300 K, 1400 K, 
1500 K, and 1600 K were obtained on the new backscat- 
tering spectrometer ^] of the NIST Center for Neutron 
Research in Gaithersburg having a sub-/ieV energy reso- 
lution and, in the set-up used, an accessible wave number 
range of q — 0.6 — 1.6 A 

The raw data reduction consists of normalization to 
a vanadium standard, correction for self absorption and 
container scattering, and interpolation to constant wave 
numbers q in order to obtain the scattering law S(q,uj). 
Further, S(q, to) was symmetrized with respect to the en- 
ergy transfer hu by means of the detailed balance factor. 
Whereas scattering from the Si and O atoms is exclu- 
sively coherent, Na scatters coherently and incoherently. 
The incoherent scattering from Na reflects itself in a flat 
background in the elastic structure factor S(q,uj = 0). 

Figure 1 displays the elastic structure factor of glassy 
and viscous sodium disilicate. Toward small q the signal 
is dominated by incoherent contributions of the sodium 
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atoms. The maximum around ~ 1.7 A -1 reflects the 
coherent scattering mainly of the disrupted tetrahedral 
Si-0 network. The intensity at q values larger than 
~1.1 A -1 is decreasing with increasing temperature and 
q which is well accounted for by the Debye-Waller fac- 
tor. In contrast, at q ~ 0.9 A -1 a pronounced prepeak 
is emerging with increasing temperature. We note that 
the signal is independent on thermal history indicating a 
process in thcrmocquilibrium as a cause for this change. 

Our S(q,u = 0) at 400 K and 800 K are consistent 
with neutron diffraction results on sodium tetrasilicate 
Na2Si40g H where at temperatures up to some 180 K 
above T g spectra display only a little change in the struc- 
ture factor. 29 Si nuclear magnetic resonance (NMR) 
measurements on sodium disilicate jl(J in the metastable 
melt up to 876 K report a continuous but small change in 
the number of O atoms surrounding a Si atom with in- 
creasing temperature. At higher temperatures motion 
of the atoms becomes too fast for conventional NMR 
techniques. From wide angle X-ray diffraction on vari- 
ous sodium and potassium silicates and sodium disilicate 
at 295 K and 1273 K in particular, previous workers |ll]] 
came to the conclusion that the structure of silicate melts 
is insensitive to temperature change. However, data qual- 
ity especially in the q range below ~ 1.2 A -1 is rather 
poor. The molecular dynamics simulations on sodium 
disilicate melts 0] do indeed exhibit a shoulder in the 
static structure factor at q — 0.94 A -1 . The molecular 
dynamics simulations Q link this prepeak to a charac- 
teristic length scale of regions where the network is dis- 
rupted and where the sodium concentration is enhanced. 

Figure 2 shows high-energy resolution backscattering 
spectra of sodium disilicate melts at g=1.6A -1 . A mea- 
surement of the glassy sample at 300 K represents the 
instrumental resolution profile R(q,u>). Above 1400 K a 
small broadening of the signal can be resolved. Around 
q ~ 1.6 A^ 1 the signal is dominated by the scattering of 
Si and O. Structural relaxation of the partially disrupted 
tetrahedral Si-0 network leads to a quasielastic broaden- 
ing in the spectra. A common feature of structural relax- 
ation in glass forming liquids is a stretching of correlation 
functions over a wider time range than expected for ex- 
ponential relaxation [ |l2|]T3]| . The line is a fit with the 
Fourier transform of the Kohlrausch stretched exponen- 
tial function convoluted with the instrumental resolution 
function: 



R(q,uj) ®{A( dte 



exp[ 



it/T q )f><])+B5(u)}. (1) 



A best fit of expression (1) with B = gives a relaxation 
time of 4 ns and a stretching exponent ft — 0.75 ±0.1. 
B = demonstrates that the dynamics of all atoms con- 
tributing to the quasielastic signal are on the same time 
scale. Scattering on considerabely faster relaxation or 
vibration (ps time scale) results in a flat background. 
At 1600 K sodium disilicate has a viscosity of 4.5 Pas 



H . The relaxation time is consistent with the relaxation 
time estimated from viscosity showing that the data in 
Fig. 2 indeed display the network dynamics. A fit with a 
stretched exponential function to the coherent intermedi- 
ate scattering function of Si and O correlations at 1.6 A -1 
as obtained by the molecular dynamics simulations [^] 
gives [3 = 0.83 in accordance with our result. 

At temperatures up to ~ 2 x T g the diffusive dynamics 
of the alkali atoms in silicate melts is appreciably faster 
than the structural relaxation G3 . Figure 3 displays the 
scattering law S(q, to) of sodium disilicate as obtained on 
the neutron time-of-flight spectrometer IN 6. Well below 
q = 0.5 A -1 the signal is dominated by the incoherent 
scattering on sodium. Diffusive motion of sodium causes 
a broad quasielastic signal. The lines in Fig. 3 are fits 
with expression (1), where the 6 function represents the 
elastic scattering contribution to the signal and B its am- 
plitude. Best fits have been obtained with a stretching 
exponent (3 within [0.7,0.8] in accordance with the MD 
simulations result of (3 — 0.77 for the incoherent interme- 
diate Na scattering function ft < 1 indicates that 
Na diffusion is not a simple activated process. Result- 
ing relaxation times (Fig. 3) show a small change with 
temperature with r ~ 29 ps at 1200 K and r ~ 13 ps at 
1600 K. 

Compared to other multicomponent glass forming liq- 
uids investigated so far the changes in the elastic struc- 
ture factor of sodium disilicate melts are unexpected. 
Viscous metallic ZrTiNiCuBe melts show a similar tem- 
perature dependence of viscosity p5| as sodium disili- 
cate when plotted versus T g /T at temperatures well be- 
low the liquidus. In addition, diffusive dynamics of var- 
ious atoms, that are smaller than the large structural 
Zr atoms, are appreciably faster than structural relax- 
ation and also decoupled from viscous flow [p^t . Toward 
higher temperatures the temperature dependence of the 
viscosities of ZrTiNiCuBe and sodium disilicate deviate 
[f§ which may very well be the result of the sodium rich 
regions in sodium disilicate: The formation of sodium 
rich regions leaves regions with a lower sodium concentra- 
tion and hence with a less disrupted Si-0 network leading 
to a higher viscosity of the entire system. 

Besides a non-exponential structural relaxation viscous 
metallic melts exhibit a fast relaxation process [|l7| which 
can be visualized as a rattling of the atoms in the tran- 
sient cages formed by their neighbors. This - so called 
- fast ft relaxation, that prepares structural relaxation 
responsible for viscous flow, has been predicted by the 
mode coupling theory (MCT) of the liquid to glass tran- 
sition |n| and has been found in many experiments and 
simulations on a variety of other glass forming systems 
p8[ . Molecular dynamics simulations on sodium disil- 
icate show that the predicted factorization property for 
Si and O correlations holds in the fast ft relaxation regime 
0. The signal in the fast ft relaxation regime, typically 
around 1 meV, is in sodium disilicate governed by the 
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incoherent contributions of the diffusing sodium atoms. 
However, since the sodium disilicate spectra do not show 
an indication for an additional fast process, its amplitude 
has to be small compared to the quasielastic line origi- 
nating from the sodium diffusion. This compares well 
to results of a neutron scattering experiment on another 
network forming melt: In liquid GeC>2 a fast MCT j3 re- 
laxation process could not be detected up to T = 2 x T g 




A recent molecular dynamics simulations study on 
sodium tetrasilicate suggests the existence of sodium 
rich channels having a distance of about 6-8 A. The pre- 
peak at 0.9 A -1 corresponds to the distance between 
these channels. A similar behavior has been found in 
the MD simulations on sodium disilicate 0. The mi- 
croscopic picture given by the molecular dynamics sim- 
ulations P,f7pC|] of regions that exhibit a higher sodium 
concentration is in agreement with our results. 

In conclusion, we have investigated structure and dy- 
namics in sodium disilicate melts with inelastic neutron 
scattering. The viscous melt well above T g exhibits 
changes in the elastic structure factor: With increasing 
temperature a pronounced prepeak emerges at ~0.9 A -1 
caused by sodium rich regions in the disrupted Si-0 net- 
work. High resolution quasielastic neutron scattering 
shows structural relaxation of the partially disrupted Si- 
O tetrahedral network to be on a nanosecond time scale 
at 1600 K. The structural relaxation is non-exponential 
with a stretching exponent f3 ~ 0.75. Inelastic neutron 
scattering reveals diffusive dynamics of the sodium ions 
on a 10 ps time scale. Even in the equilibrium melt the 
Sodium ions diffuse in a relative immobile Si-0 matrix. 
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FIG. 1. Elastic structure factor S(q,co = 0) of sodium dis- 
ilicate: Toward small q the signal is dominated by the in- 
coherent contributions of the sodium atoms. The maximum 
around ~ 1.7 A -1 reflects the scattering of the partially dis- 
rupted tetrahedral SiCU network. With increasing tempera- 
ture a pronounced shoulder is emerging at g~0.9A _1 . The 
decrease in intensity at q values larger than ~ 1.1 A -1 in turn 
mainly reflects the Debye-Waller factor. 
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FIG. 2. High-energy resolution backscattering spectra of 
a sodium disilicate melt. Around q ~ 1.6 A -1 the signal is 
dominated by the coherent scattering of Si and O. At 1600 K 
structural relaxation of the Si-0 network is on a ns-time scale 
causing a broadening of the quasielastic line. The line repre- 
sents a fit (1) giving a relaxation time of ~ 4 ns and a stretch- 
ing exponent /3 = 0.75 ±0.1. 
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FIG. 3. Scattering law S(q,u>) at q = 0.5 A -1 measured 
on the neutron time-of-fiight spectrometer IN 6. Diffusion of 
sodium ions results in a broad quasielastic signal that is de- 
scribed by (1) with a stretching exponent (3 within [0.7,0.8] 
and relaxation times r at some 29 ps at 1200 K and some 13 ps 
at 1600 K. 
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